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ABSTRACT 

Active galactic nuclei (AGN) have Fe ii emission from the broad line region (BLR) that differs greatly 
in strength from object to object. We examine the role of the total and gas-phase iron abundance 
in determining Fe ii strength. Using AGN spectra from the Sloan Digital Sky Survey (SDSS) in the 
redshift range of 0.2 < z < 0.35, we measure the Fe/Ne abundance of the narrow line region (NLR) 
using the [Fe vii]/[Ne v] line intensity ratio. We find no significant difference in the abundance of Fe 
relative to Ne in the NLR as a function of Fe il/H/3. However, the [N ii]/[S ii] ratio increases a by 
a factor of 2 with increasing Fe ii strength. This indicates a trend in N/S abundance ratio, and by 
implication in the overall metallicity of the NLR gas, with increasing Fe ii strength. 

We propose that the wide range of Fe ii strength in AGN largely results from the selective depletion of 
Fe into grains in the low ionization portion of the BLR. Photoionization models show that the strength 
of the optical Fe ii lines varies almost linearly with gas-phase Fe abundance, while the ultraviolet Fe ii 
strength varies more weakly. Interstellar depletions of Fe can be as large as two orders of magnitude, 
sufficient to explain the wide range of optical Fe ii strength in AGN. This picture is consistent with 
the similarity of the BLR radius to the dust sublimation radius and with indications of Fe ii emitting 
gas flowing inwards from the dusty torus. 
Subject headings: galaxies: active — quasars: general 



1. INTRODUCTION 

The broad emission-line spectrum of quasars often in- 
cludes strong Fe ii in the optical and ultraviolet. The dif- 
ference between the weakest and strongest optical Fe ii 
emission exceeds a factor of 10, measured as equivalent 
width (EW) or as Fe ii/H/3 line ratio. For recent discus- 
sions of Fe II intensities in AGN, an d refe rences to ear- 
lier w ork, see iKovacevic et all ()2010D and iFerland et al.l 
()2009t) . The Fe li strength anti-correlates with the 
strength of the narrow [O ill] emission line. This trend 
along with several associated correlations defines the so- 
called "Eigenvector 1" (EVl), which characterizes some 
of the m ost conspicuous differenc es among the properties 
of AGN (jBoroson fc Green|[l992l hereinafter BG92). The 
quest for physic al drivers of EVl has inspired a numbe r 
of studie s (e.g IWills et al.l 119991: IMa rziani e t al.l 120031) . 
IBorosonI (|2002| ) and INetzer fc Trakht enbrot ( 2007f ) find 
that Fe ii increases with Eddington ratio L/LEd, as orig- 
inally suggested by BG92. However, the physics under- 
lying this correlation remains unclear. There is even de- 
bate as to whether the Fe ii emission is entirely pow- 
ered by the ionizing continuum of the central source, 
or comes in some measure from a mechanically heated 
region (IWills . Netzer. & Wills 1 985; CoUin .& Jolv 200Ci; 
ISigut fc Pradhani 120031: iBruhweiler fc Verneii [20081) . In 
this situation, any observational clues to the nature of 
EVl and the great range of Fe ii strength are of value. 

It is generally assumed that the strength of Fe II emis- 
sion is driven by physical conditions within the BLR, 
such as ionizing continuum, BLR density and geometry. 
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column density, and turbulent velocity. However, a high 
abundance of Fe has been disc ussed to help produce 
the strong est Fe ii observed (IWills. Netzer. fc Willp 
1985t ICoUin-Souffrin. Hameurv. fc Jolvl [TgM 



Hamann fc Ferlancilll993D . The utilitv of Fe ii to assess 



the Fe abundance in high redshift QSOs has received 
considerab le interest, in the context of galactic chemical 
evolution (iHamann fc Ferlandl[l9 9^; Verner et al.l JlOOl 
iBaldwin et al.l 120041 : INetzer fc Tra khtenbrot 2003, and 
references therein) . Here we assess the importance of the 
abundance of Fe relative to the a-elements, the overall 
metallicity of the nuclear gas, and the depletion of Fe 
into grains for the strength of Fe li emission in QSOs. 

2. QSO SAMPLE AND MEASUREMENTS 

We investigated the influence of differing abundances 
in quasars on their optical Fe il emission strength by 
studying the optical emission-line properties of a sample 
of QSOs from the Sloan Digital Sky Survey (SDSS) B 
It is difficult to measure abundances within the BLR di- 
rectly, because it is a region of high density, line width, 
and line optical depth. On the assumption that the abun- 
dance in the NLR and BLR is the same for a given object 
(see below), we used the narrow emission-line spectrum 
to assess abundances in the NLR. We considered two 
key line ratios: (1) The [Fe vii]/[Ne v] intensity ratio 
gives a measure the Fe/Ne abundance ratio. (2) The 
[N ii]/[S ii] ratio gives a measure of the N/S abundance 
ratio, which is in turn a secondary indicator of the over- 
all metallicity of the gas. Our sample consists of 1571 
quasars from SDSS Data Release 7 (DR7). These ob- 
je cts were selected i n the manner of the "HOB" sample 
of iSalviander et all (j2007f ). with the additional require- 
ment of a signal-to- noise (S/N) ratio greater than 10 in 

^ The SDSS website is |http: / /www.sdss.org| 
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the continuum at A5100 rest wavelength. The flux and 
equivalent width (EW) of the broad H/3 line and the opti- 
cal Fc II blends were measur ed with the aid of a spec trum 
fitting program described bv lSalviander et all (j2007l ). us- 
ing a template fitting procedure to establish the flux in 
Fe II relative to the local continuum. We characterized 
the Fe ii emission strength using the flux ratio of the 
Fe II 4570 A blend to broad H/3, following BG92. We 
used quasars at redshifts 0.2 < z < 0.35 to ensure cov- 
erage of both [Ne v] A3425 and [Fe vii] A6087. In order 
to bring out the weak [Fe vii] line, we binned the 1571 
objects by Fe li strength, and made five composite spec- 
tra of "very weak," "weak," "medium," "strong," and 
"very strong" Fe il emission. These composites had 312 
to 315 objects, within bins bounded by Fe ii/H/3 values 
of 0, 0.215, 0.372, 0.524, 0.708, and 1.82, respectively. 
Individual spectra were corrected for Galactic redden- 
ing using the extinction values Ag given by the SDSS 
pipeline, and normalized to a mean flux density F\ of 
unity using all wavelength points in a particular spec- 
trum. The individual spectra in each group were shifted 
in wavelength to the rest frame and re-binned to a com- 
mon wavelength grid at a spacing of 1.41 A. The adopted 
composite spectrum was a mean of the rebinned F\ for 
all contributing spectra at a given wavelength. The com- 
posite spectra are shown in Figures [T] and [5] The region 
of the [Fe vii] and [Ne v] line is shown in Figure [2] For a 
dis cussion of issues involving composite spectra of QSOs, 
see lVanden Berk et al.l (|2001| ). and references therein. 

From these composites, we measured the emission-line 
fluxes of a number of lines, including [O ill] A5007, 
[O II] A3727, [Ne iii] A3869, [Ne v] A3425, 
[S ii] AA6716, 6730, and [Fe vii] A6087. The results are 
given in Table 1 and Figure [31 where the Fe ii/H/? val- 
ues are averages of the values for the individual spectra 
that compose each composite. Most lines were measured 
using a Gaussian fit with the IRAF task SPLOT0. The 
broad Balmer emission lines (Ha, H/3) were measured 
using a Lorentzian profile. Estimated uncertainties are 
10% for the stronger lines, including continuum place- 
ment and faithfulness of the fit. For [Fe vii] and [Ne v], 
the uncertainty is as much as 20%, based on noise, con- 
tinuum uncertainty, and the presence of a strong blue 
wing on both lines that we excluded from the fit. The 
[N ii] A6583 and A6548 lines were fairly weak bumps on 
the wings of the broad Ha line and relatively difficult to 
measure. Therefore we measured the [N ii] intensity by 
subtracting from the Ha - [N il] blend a doublet with 
the theoretical 3.0-to-l intensity ratio, each line having 
a Gaussian profile with a central wavelength and width 
based on the redshift and line width of [S ii] A6716. The 
intensity of the doublet was adjusted so that the Ha line 
had a smooth profile with no visible residual intensity or 
over-subtraction of [N ii] . Error bars were estimated by 
determining [N ii] intensities giving a slight under- or 
over-subtraction as judged by eye. This gave an uncer- 
tainty of about ±12% for each composite. 

Our results will be discussed in terms of trends of ob- 
served line ratios. For actual ionic abundances, collision 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



TABLE 1 

Emission Line Ratios for Composite Spectra 



Emission 






Intensity Ratio 








Lines 


vlow 


low 


med 


high 


vhigh 


N II 


/ 


S II] 


1.21 


1.25 


1.46 


1.62 


2.15 


N n 


/ 


O II] 


0.96 


1.25 


1.22 


1.55 


2.29 


O II 


/ 


O ml 


0.24 


0.23 


0.24 


0.25 


0.26 


[S ii]/[0 II] 


0.80 


0.95 


0.83 


0.95 


1.07 


Fc vii]/[Ne v] 


0.33 


0.32 


0.37 


0.35 


0.32 


Ne v]/[Ne in] 


0.97 


1.10 


1.15 


1.61 


1.81 


Ho/H/3 


3.64 


3.39 


4.14 


3.52 


3.04 


"S II] (6720/4072) 


9.2: 


10.3: 


9.7: 


9.9: 


13.5: 


"S ii] (6716/6731) 


1.19 


1.05 


1.07 


1.10 


1.08 


Fe ii/H/3 


0.11 


0.29 


0.44 


0.61 


0.91 



Note. — Intensity ratio for emission lines measured from com- 
posite spectra binned by broad Fe II strength. Values refer to the 
narrow emission lines except for the Balmer lines and Fe II. In- 
tensities include both lines for the [O ll] and [S ll] doublets but 
only the stronger line for [N ll], [O III], [Nc III], and [Ne v]. Colon 
indicates large uncertainty. Sec text for discussion. 

strengths from Bcrrington ct al. (2000) and Osterbrock 
& Fcrland (2006) lead to the relation 

n{Fe+^)/n{Ne+'^) = 0.91 /(A6087)//(A3425) (1) 

for an assumed T=15,000 K based on photoionization 
models. Likewise, for the [N ii]/[S ii] ratio we have 

n{N+)/n{S+) = 6.0 /(AA6584,6548)//(AA6716, 6730) 

(2) 

for an assumed T=12,000 K. These expressions give 
the ionic abundance ratios, if collisional de-excitation is 
unimportant. 

3. RESULTS FOR THE NARROW LINE REGION 
3.1. Iron 

iNussbaumer fc Osterbrock ()197Gf ) suggested that 
[Fe vii]/[Ne v] should be a good measure of the Fe/Ne 
ratio, based on the similarity of the ionization potentials. 
There is no significant trend in [Fe vii]/[Ne v] in our 
composite spectra. The Fe ii/H/3 ratio varies by a factor 
of 8 from the "very low" to "very high" composite. 
Taking neon to represent the a-elements in general, we 
conclude that differences in Fe abundance, relative to 
the a-elements is not a significant cause of the observed 
range of Fe ii emission strength, from a statistical point 
of view. In particular, overabundances of Fe/0 and 
Fe/Mg of a factor of 2 to 10, as motivated by chemical 
evolution models and attempts t o fit the Fe ii/JVIg ii ra- 
tio with photoionization model s (j Wills. Netzer. fc Willsl 
Il985t iHamann fc Ferlandlll993l ). appear to be ruled out. 

This conclusion is based on the assumption that the 
[Fe vii]/[Ne v] ratio is a faithful measure of Fe/Ne, and 
that Fe is not significantly depleted into g rains i n the 
[Fe vii] zone of th e NLR. iFerguson et al.l ()1997aD and 
iNagao et al.l (|2003[ ) conclude that refractory elements are 
not depleted in the coronal line region of the NLR, but 
they also find that [Fe vii] and [Ne v] do not come 
from the same place in "locally optimally emitting cloud" 
(LOG) models of the NLR. Here we assume that any 
ionization correction for Fe+^/Ne"'"^ in the NLR does 
not change systematically with the broad line Fe ii in- 
tensity among our composites. We further assume that 
the abundances in the BLR are similar to those in the 
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NLR. However, it is possible that intense star formation 
in the nucleus may give chemical e nrichment on a spa- 
tial scale smaller than the NLR (see lHamann "erall[2007l: 
iHamann fc Ferlandlll999l . and references therein). 

3.2. Ionization 

The composite spectra were originally constructed to 
assess the [Fe vii] strength in the NLR, but they also 
afford an opportunity to examine other narrow line ra- 
tios for any systematic dependence on Fe ii strength. 
One issue is the level of ionization of the gas. Tabic 1 
and Figure |3] show that the [O ii] / [O ill] ratio is closely 
similar for the various Fe ii bins. There is a significant 
increase in [Ne v]/[Ne ill] with increasing Fe ii strength, 
reflecting a ^ 25% decrease in the [Ne v] equivalent width 
across the bins and a much larger decrease in the EW of 
[Ne III]. The narrow He ii A4686 line to [O ill] A5007 
intensity ratio (not given in Table 1) shows a similar, 
though less continuous, increase from 0.06 to 0.13 across 
the five bins; this reflects a substantially constant equiv- 
alent width of narrow He ii together with a systematic 
decrease in the EW of A5007 from 25 to 12 A across the 
bins. The anticorrelation of Fe li and [O ill] is well known 
(BG92). Although the constancy of [Fe Vll]/[Ne v] could 
result from offsetting ionization and abundance trends, 
a straightforward interpretation is that the relative size 
of the zone containing O"*"^ and Ne"*"^ decreases with in- 
creasing Fe II without a major effect on the Fe+^/Ne+"' 
ionization correction. 

3.3. Reddening 

The [S ii]/[0 il] ratio shows an increase of about 30% 
with increasing Fe ii strength. This may be an indica- 
tion of a modest increase in reddening of the NLR with 
increasing Fe ii. The [S ii] /(A6720)//(A4072) ratio is 
quite uncertain but is consistent with a reddening trend 
of this magnitude. This trend is not evident in the broad 
/(Hq!)/I(H/3) intensity ratio. 

3.4. Nitrogen 

The [N ii]/[S ii] ratio shows a systematic increase by 
a factor of 2 from the very low to the very high Fe ii 
bins. This ratio is insensitive to reddening and electron 
temperature. The N"*" and S+ ions have similar ionization 
potentials and occupy similar zones of the nebular ioniza- 
tion structure. This and the constancy of [O ii]/[0 ill] 
suggests that the trend in [N ii]/[S ii] is not a result 
of ionization of the NLR. The trend of increasing [N ii] 
strength is evident in the [N ll]/[0 ll] and [N ii]/[0 ill] 
ratios as well. These results imply a real trend in the 
N/S and N/0 chemical abundance ratios with increas- 
ing Fe II strength, amounting to a factor of 2 from the 
"very low" to the "very high" Fe ii strengths. Nitro- 
gen is largely a secondary nucleosynthetic product, so 
th at N/0 increases wi th 0/H. The H ii region results 
of Ivan Zee et~all (|l998f ) show N/0 increasing almost lin- 
early with O /H above 12-t-log O /H = 8.5. Assuming that 
0/H is in this range in the AGN studied here, then the 
trend of N/0 with Fe ii implies an increase of a factor 
of 2 in 0/H. Although the chemica l evolution of AGN 
host galaxies may be complicated pTamann fc FerlanJ 
[T99llHamann et al.|[200llNetzer fc Trakhtenbrotll2007D . 
it may be reasonable to assume that Fe/H varies roughly 



with 0/H for the modest redshifts considered here. In 
this case, the Fe/H abundance ratio may contribute 
roughly a factor of 2 to the range in Fe ii strength from 
our "very low" to "very high" bins. This conclusion 
is con sistent with the work of iNetzer fc Trakhtenbrod 
(|2007f l. who find an increase of Fe ii strength with Ed- 
dington ratio L/Lsd in a large sample of SDSS quasars. 
Combining this with in dications of an increasing N/C 
with L/^Ed in the BLR (jScott et al.ll2004l: IHamann et al.l 
I2OOI ). they argue that the overall metallicity in QSOs in- 
creases with Fe II strength. However, our results indicate 
that this makes only a modest contribution to the full 
range of optical Fe 11 strength in AGN. Note, however, 
that our QSO sample is at low redshift, whereas much of 
the interest in iron abundances in QSOs has focused on 
high redshifts. 

4. Fe II STRENGTH AND X-RAY HEATING 

It is widely assumed that the Fe il emission is largely 
powered by the soft X-ray portion of the ionizing spec- 
trum, given that photoionization is indeed the primary 
excitation mechanism. The harder photons in the ioniz- 
ing continuum create an extensive warm, partially ion- 
ized zone where Fe 11 and other low ions are subject to 
coUisional excitati on and in some c ases o t her excitation 
mechanisms (see iKwan fc KroTig Il98ll: iFerland et afl 
120091 and references therein) . Does the relative strength 
of the X-ray continuum drive the range of Fe 11 strength 
observed in AGN? The following considerations suggest 
that it does not: 

1) A number of empirical studies have examined corre- 
lations of Fe II stren gth with the X-ray slope or with the 
X-ray/optical ratio (jSulentic. Marziani. fc Dultzin|[2000l 
and references therein). Lawrence et al. (1997) stud- 
ied a sample of AGN with extreme values of R-pc 11 = 
/(A4570)//(H/3). The strong Fe 11 emitters have X-ray 
properties, including ofox, similar to the weak Fc 11 ob- 
jects. In particular, the prototype strong Fe 11 object, 
I Zw 1, has Q!ox = —1-4, an entirely typical v alue. Com- 
bining their data with the complete sample of lLaor et al.l 
(|1997f l. Lawrence et al. found little correlation between 
aox and i?Fe 11, and only an ambiguous correlation with 
Qfx- Indeed, Figure 3 of Lawrence et al. shows a slight 
trend in the sense of strong optical Fe 11 for weak X- 
ray luminosity (steep ctox)- They did find a significant 
correlation with the X-ray to IR slope ctix in the sense 
of stronger infrared for stronger Fe 11. Using co mposite 
spectr a for X-ray bright and X-ray faint QSOs, iGreenI 
(fl99l found that UV Fe II was stronger whereas opti- 
cal Fe II was weaker for X-ray bright objects. However, 
the differences were small compared to the full range 
of Fe II strength among individual AGN. iLeighlv et al.l 
(|2007l) find a weak X-ray continuum but strong Fc 11 
emission in PHL 1811, and discuss other similar exam- 
ples. None of these results gives support for the idea 
that stronger X -ray continuum drives str onger Fe II for 
AGN as a class. IFerland fc PerssonI ()1989[ ) reached a sim- 
ilar conclusion regarding the Ca 11 emission from AGN, 
which also comes from the partially ionized zone. 

2) In order to explore the expected response of the 
Fe II emission to differences in X-ray luminosity (and 
other parameters) we have computed a set of models of 
the BLR using version 07.02.00 of the photoioni zation 
code Cloudy, most recently described by .Ferland et al.l 
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TABLE 2 

Photoionization Model Results 



Model 


1 


2 


3 


4 




-1.4 


-2.0 


-1.4 


-1.4 


Fe Depletion 


0.0 


0.0 


-1.5 


-1.5 


Other Depletion 


0.0 


0.0 


0.0 


-1.5 












ioUU 


O.z 


1 .0 




i.oy 


2500 


17.8 


21.0 


5.1 


3.9 


OUUU 


U.OO 


U.OO 


n non 
U.Uzy 


U.UZ4 


7000 


0.069 


0.081 


0.0059 


0.0045 


10000 


0.42 


0.44 


0.030 


0.024 


Mg 11 


3.4 


3.2 


2.8 


0.20 


Si III] 


1.30 


1.46 


1.19 


0.037 


Si IV 


1.56 


1.75 


1.42 


0.104 


Hq 


4.6 


4.6 


4.5 


4.6 


Lya 


28.1 


32.4 


25.8 


21.6 


C III] 


1.69 


1.78 


1.54 


1.45 


C IV 


13.8 


13.2 


12.6 


11.1 



Note. — Results of Cloudy photoionization models. 
Logarithmic depletion relative to hydrogen is given for 
iron and for other refractory elements. See the text for 
discussion. 



1998) 



As a reference model, following iFerland et alj 
2009D, we used solar abundances, an ionizing flux ^ = 

The 
and 



10^^ cm ^ s ^, and a gas density N = 10^° cm 



ing an ionization parameter U = (j)/Nc = 10 
internal turbulent velocity was Mturb = 100 km s~'^ 



the stopping column density was 10 cm"'', in order 
to include an extensive partially ionized zone. The ion- 
izing continuum was the sum of (1) a UV component 

(X v-°-^cxp{-hiy/kTcut) with T^nt = lO^ "^ K to rep- 
resent the Big Blue Bump, and (2) an X-ray power law 
Li, oc h'~^. The continuum had a low frequency expo- 
nential cutoff below 0.01 Ryd. The models used the full 
treatment of the Fe ii ion (371 levels) and were iterated 
to convergence of the diffuse radiation field. The ratio of 
the X-ray to the UV component is controlled by the pa- 
rameter ttox- For the reference model (M odel 1), we used 
Oox = —1-4, a typical observed value (jLawrence et alJ 
Il997| ). This model has not been optimized to fit a typ- 
ical AGN emission-line spectrum, but simply serves as 
a reference point for exploring the effect on Fe ii of 
changing various model parameters. A simultaneous fit 
to AGN emission-line spectra requires a combination of 
photoionized clouds with a range of physical conditions 
ijBaldwin et"an [19951 ) . 

Table 2 gives line intensities from the Cloudy mod- 
els. The Ha, C ill], C iv, and Mg ii intensities are rea- 
sonable. Lya is stronger than observed relative to H/3, 
a familiar problem w ith photoionization models (e.g., 
IKwan fc Kroliklll981[) . The Cloudy output gives Fe ii 
intensities summed over broad wavelength bands at 1000 

- 2000, 2000 - 3000, 4000 - 6000, 6000 - 7800, and 7800 

- 30000 A. Note that the A2500 band includes the broad 
UV Fe II bump, and the A5000 band contains the promi- 
nent optical Fe li blends at A4570 and A5250. Model 1 
under-predicts the Fe ii intensity seen in strong Fe II ob- 
jects, a common problem with photoionization models as 
mentioned above. 



In Model 2, we used aox = —2.0 to explore the ef- 
fect of weaker X-rays while remai ning within the span o f 
frequently observed values of aox (jLawrence "eFaII[l997l) . 
The Fc II intensity (relative to H/3) actually increased 
slightly in the X-ray weak model. We also computed 
an alternative pair of models with </> = 10^^ cm~^ s^^, 
N = 10^ cm-3, Tcut = 10^-^ K, and no turbulence. In 
this case, the Fe ii A2500 blend decreased by 0.16 dex and 
the A5000 band by 0.20 dex in going from aox = —1.4 to 
—2.0. Even in this case, the degree of change is insuffi- 
cient to give the full observed range of Fe ii strength. An 
extreme systematic variation of aox with Fe ii would be 
required, for which the observations give little support. 

5. Fe II STRENGTH AND IRON DEPLETION INTO GRAINS 

The above results indicate that changes in the total ele- 
mental abundance of Fe and the X-ray luminosity do not 
cause the wide observed range of Fe ii emission strength 
among AGN. We propose instead that Fe ii strength in 
AGN largely results from differing degrees of depletion 
of the gas-phase abundance of iron into grains in the 
relevant portions of the BLR. Gas-phase depletions of 
refractory elements in the interstellar medium can be se- 
vere. Iron is depleted by up to 2 orders of magnitude in 
the interstellar medium and in ionized nebulae (see dis- 
cussion below). Such a degree of depletion, if present in 
AGN with weak Fe II but not those with strong Fe ii, can 
account for much of the observed range of Fe ii strength 
in AGN. 

There has been considerable discussio n of refractory 
element depletions in the NLRs of AGN (iFerguson et al.l 
ll997bD . IGaskell. Shields, fc Wam^ ((19811) considered 
the question of refractory element depletions in the BLR 
and concluded that depletions of Si, Mg, and Fe as severe 
as those in the ISM did not occur. This conclusion was 
based mostly on the Mg ii. Si ill, and Si iv lines, and 
the intensity of Fe ii in strong Fe li objects. Here we 
suggest that the degree of depletion of iron varies from 
object to object, and can be severe in objects with weak 
Fe II. In a study of I^Q i] and Can emission from AGN, 
iMatsuoka. Kawara. fc Ovabul ()2008f ) note that depiction 
of Ca II might help to reconcile predicted and observed 
intensities and mention the possibility that deple tions 
may affect the Fe ii lines. IFerland fc PerssonI ()1989() note 
that dust mixed with the BLR gas could help to explain 
low observed values of the Ca ii H and K emission lines, 
relative to the infrared triplet. As discussed below, the 
radius of the BLR is interestingly close to the dust sub- 
limation radius; and for typical parameters, grains may 
exist in the partially ionized zone of the BLR clouds but 
not in the highly ionized surface layers. 

The present proposal does not address the long- 
standing question of how to e xplain the large Fe ii 
strength observed in many AGN ()Wills. Netzer. fc Willj 
ll985HBruhweiler fc Verneiil2008l and references therein). 
Rather, it serves to separate the question of Fe II strength 
into two parts: (1) why is Fe ii so strong in some ob- 
jects, and (2) why is it so weak in others? We do not 
attempt to resolve the first question in this paper, ex- 
cept to note that under the present proposal, the Fe ii 
physical excitation mechanism is freed from the require- 
ment of explaining by itself the wide range of observed 
Fe II strengths. 
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5.1. Photoionization Models 

In order to explore the dependence of Fe ii strength 
on gas-phase Fe abundance we have computed two ad- 
ditional Cloudy models (see Table 2). Model 3 is oth- 
erwise identical to Model 1 but has a depletion of iron 
by 1.5 dex. The A2500 feature decreased in intensity 
by 0.66 dex, and the A5000 feature by 1.36 dex, corre- 
sponding to (Fe/H)°'^"' and (Fe/H)°-^^, respectively. In 
the N = 10^ cm~^ model described above, a 1.5 dex de- 
pletion of Fe abundance reduced the A2500 and A5000 
features by 0.82 dex and 1.37 dex, respectively. Thus, 
the optical Fe ii intensity varied almost linearly with the 
gas-phase iron abundance, and the UV Fe ii less strongly. 
These results a,re consistent with p revio us photoioniza- 
tion s tudies. iVerner et al.l ()2003l) and [Baldwin et al.l 
(|2004[ ) found that the Fe il intensity varies approximately 
as (Fe/H)° * for the optical lines and as (Fe/H)"-'' for 
the UV bump. Some of these models assume differ- 
ent physical p aram eters for the BLR from ours, and the 
IVerner et all (|2003| ) models use an 830 level Fe"*" model 
at om as opposed to 37 1 levels in our models and those 
of [Baldwin erall ([200l . 

Other refractory elements may be depleted where iron 
is. In Model 4, the abundances of Al, Mg, Si, Ca, and 
Fe were all depleted by 1.5 dex. Table 2 shows that the 
multiplets Mg ii A2800, Si ill A1892, and Si iv A1400 
decreased roughly in proportion to the gas-phase abun- 
dance. The behavior of the Fe ii multiplets in Model 4 
is similar to the case in which only Fe was depleted. 

Different refractory elements have different depletions 
in the ISM, and different grain compositions may be more 
or less easily destroyed in the AGN environment. This 
may lead to useful di agnostics for grain destructio n and 
gas-phase depletions. iDelgado Inglada et al.l (120091 ) sum- 
marize published depletion factors in planetary nebulae: 
1/6 to 1/300 for Ca, 1/2 to 1/350 for Al, 1/3 to 1/300 for 
Fe, near solar to 1/10 for Mg, and near solar to 1/20 for 
Si. Silicon appears to be depleted by a le sser factor than 
iron in H ii regions. IGarnett et all (|l995[) found Si deple- 
tions of only -0.1 to -0.6 dex in extragalactic H II regions, 
significantly less than in dense interstellar clouds. If Si is 
more easily restored to the gas phase in H II regions, the 
same may be true in AGN, allowing objects with weak 
Fe II to have normal intensities of Si ill and Si IV. 

5.2. Geometrical Considerations 

What might cause widely differing degrees of depletion 
of iron, and possibly other refractory elements, among 
AGN? From considerations of BLR covering factor, Ly- 
man continuum absorption, line widths, and reverbera- 
tion mapping, GaskcU (2009;) argues that the BLR is the 
inward extension of the dusty torus, with an inflow veloc- 
ity a substantial fraction of the orbital velocity. Such a 
picture lends itself to the idea of refractory element deple- 
tions in the low ionization zone. As material flows inward 
through the dusty torus, its equilibrium temperature 
rises as it experiences a stronger radiation field, eventu- 
ally reaching the sublimation temperature of ^ 1500 K. 
The exact su blimation radius depe nds on grain size and 
composition ()Laor fc Draind [19931 ). but there will be a 
point where the refractory elements are substantially re- 
stored to t he gaseous pha se. This sublimation radius 

IS given bv iLaor fc Draind (|1993D as Rd = (0.2 pc)L46, 



where L^q is the AGN luminosity in units of lO"'^ erg s~^. 
If this happens before the material reaches the main re- 
gion of low ionization line emission, then strong emission 
in Fe ii and other lines of refractory elements will occur. 
On the other hand, if the dust survives through most of 
the low ionization zone, then these emission lines will be 
weak, because of the lack of emitting ions and because 
of attenuation of the ionizing radiation and line emission 
b y the dust itself . 

INetzer fc Laorl (jl993D suggested that the survival of 
dust outside the sublimation radius leads to extinction 
of the ionizing continuum and suppression of line emis- 
sion, setting a natural limit to radius of the BLR (see also 
lLaorll200"7l ). Infrared variability studies give a radius of 
the d usty torus just outside the BLR (iSuganuma et al.l 
|2006( ). supporting this picture. Here we suggest that 
the relationship between the sublimation radius and the 
BLR outer boundary differs from object to object. Crit- 
ical to this picture is the actual radius of Fe ii emission 
relative to the dust sublimation radius. Reverberation 
studies and line widths indicate that the lower ionization 
lines often come from larger radii in the BLR (jPetersonI 
[1997: Sulentic, Marziani. & Dultzin 2000). This may re- 
sult from some combination of actual ionization strati- 
fication and line emissivity effect s as considered in the 
LOG model (jBaldwin et al.|[T995[ ). There have been few 
successful reverberation studies of the Fe ii emitti ng ra- 
dius i n AGN. For the Seyfert galaxy NGC 5548, [Mao3 
(jl993[) found a similar time lag for the UV Fe ii lines 
as for Lya. In a reverb eration study of Ark 120, which 
has strong optical Fe li, iKuehn et al.l (|2008[ ) found that 
the lag for the optical Fe ii lines was ill-defined but may 
be around 300 days, larger than for H/3. They estimate 
a dust sublimation radius in this object of ^ 460 light 
days, and conclude that within the uncertainties "it is 
plausible that the optical Fe ii emission is produced at 
or just inside the dust sublima tion radiu s ." In a study of 
hue profiles of AGN in SDSS, iHu et al.l (l200l conclude 
that the Fe ii emission comes from an infiowing zone in 
the outer part of the BLR. Taken together, these studies 
are consistent with Fe ii emission from gas entering the 
BLR from the dusty torus. 

A refinement of this picture takes account of the atten- 
uation of the ionizing radiation field with depth within 
the ionized cloud or ionized surface layer of the disk 
or torus. Reverberation measurements of the BLR ra- 
dius (often involving the II/3 line) typically give a ra- 
dius smaller than the dust sublimation radius by a fac- 
tor of order 2, bas ed on radii and luminosities given by 
iBentz et al.l ()2009D . Thus, refractory grains should not 
survive at the illuminated face of a BLR cloud. However, 
the radiation field diminishes with depth in the cloud, 
allowing dust to survive at large depths. As the dusty 
gas fiows from deep in the torus toward the irradiated 
surface, the ambient radiation field intensifies and the 
grain temperature rises, reaching the sublimation point 
at some depth in the ionization structure. Grain equi- 
librium temperatures calculated by the Cloudy program 
support this picture. In the reference model (Model 1) 
described above, the temperature of silicate grains with 
a radius of 0.094 microns is 680 K at the maximum col- 
umn depth of 10^'^ cm~^. This allows grains to survive 
at this depth, even though the incident fiux at the cloud 
face would easily evaporate refractory grains. The grain 
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temperature rises with decreasing depth, reaching the 
subhmation temperature of 1500 K at a column den- 
sity of lO^i-^'^ cm~^. This is well below the "Stromgren 
depth" of lO^^-^'^ cm~^, where hydrogren is 50% ionized. 
The grain temperature as calculated by Cloudy would 
be 2000 K at this latter depth, an d 3100 K at th e cloud 
surface. Reference to Figure 3 of lFerland et"all ([20091) 
shows that in their model (similar to our Model 1), most 
of the optical Fe ii emission occurs below the sublima- 
tion depth of 10^^-^^ cm~^, whereas most of the ultra- 
violet Fe II emission occurs at shallower depths. Thus, 
evaporation of the grains occurs at a depth giving se- 
vere reduction of the optical but not the ultraviolet Fe ii 
emission. Dust will also not affect the emission from the 
highly ionized gas above the Stromgren depth, such as 
Si III] and Si iv; and weak optical Fe ii could accompany 
a normal ratio of the ultraviolet Fe ii to Mg ii lines. The 
exact sublimation depth depends on grain size and com- 
position, but the qualitative pattern remains for a range 
of compositions and sizes included in the Cloudy out- 
put. For different AGN with different parameters, the 
sublimation point wi ll occur at different depths in the 
ionization structure. iGaskell ct al. (2007) have also dis- 
cussed the affect of attenuation of the AGN continuum 
in the BLR on the radius of dust sublimation. 

An alternative geometry for the BLR involves 
a radiatively driven w ind from the accretion disk 
(jMurrav fc ChiangI Il998| ). In this picture, one issue is 
whether grains can survive in the inflowing material to 
the radius where it is expelled in th e wind. Standard 
accretion disk physics ()Petersonlll997[) gives for the disk 
effective temperature due to internal viscous dissipation 

T^ff = (10=^-12 K){M/M^f/Hl-''\%1,, (3) 

where M / Me is the accretion rate relative to that giv- 
ing the Eddington luminosity, Mg is the black hole mass 
in units of 10 M©, and W3000 is the orbital velocity at 
the radius of interest in units of 3000 km s~^. Thus, the 
sublimation temperature for refractory grains is reached 
at orbital velocities appropriate for the broad emission 
lines. This may be consistent with differing degrees of 
grain evaporation in different objects. (The mid-plane 
temperature of the disk will be higher.) The above ex- 
pression qualitatively agrees with the observed trend of 
stronger optical Fe 11 with increasing Eddington ratio, 
but the dependence on black hole mass may be prob- 
lematic. Moreover, it does not give such a natural L^'^ 
dependence for the BLR radius as do es the dust subli- 
mation model of lNetzer &: LaoJ (|1993[ ). Note that in the 
disk-wind picture, the energy for the line emission still 
comes from photoionization by the central continuum; lo- 
cally produced energy is insufficient at this shallow depth 
in the gravitational potential. Sublimation at an inter- 
mediate depth in the ionization structure, as discussed 
above, could also occur in the disk-wind model. 

5.3. Turbulent Velocity and Column Density 

Local turbulence substantially affects the Fe 11 spec- 
trum in photoionization models by facilitating continuum 
and line-line fluorescence. Increasing the turbulence can 
increase the Fe 11 strength and give better agreement be- 
tween t he predicted shape of the Fe il blends a nd obser- 
vation (jBaldwin et al..,2004t iBruhweiler fc Verner, ,20081 



and references therein) . IBruhweiler fc Vernerl ()2008D find 
a factor of 2 increase in the UV Fe li strength rel- 
ative to Mg II as the turbulence increases from 5 to 
50 km s~^ (their Table 2). When we decreased Uturb from 
100 km s"^ to zero in Model 1, the A2500 Fe 11 band de- 
creased a factor 10, and the A5000 band decreased by a 
factor of 2. If these results are representative, they sug- 
gest that differing amounts of turbulence could make a 
substantial contribution to the range of strength of Fe 11 
in AGN. On the other hand, turbulence still appears to 
be inadequate to give the full observed range in optical 
Fe II. Moreover, [Baldwin et al.l ()2004[) argue that sub- 
stantial turbulence is required to fit the detailed shape 
of the UV Fe 11 feature, so that there may be limited 
freedom to vary the turbulence. A useful test may be to 
examine the observed shape of the UV and optical Fe 11 
blends as a function of optical Fe 11 strength, for com- 
parison with photoionization models that vary either the 
turbulence or the gas-phase abundance of Fe. 

The column density of the emitting cl ouds also has 
an im portant effect on the Fe 11 emission. iFerland et al.l 
(|2009!) illustrate the increase in Fe ii/H/? with increasing 
column density. They find that the minimum column 
density is ~ 10^'^ cm~^ for gravity to overpower radiation 
pressu re and allow infall of clouds as found by iHu et ahl 
(|2008l ). The UV Fe 11 lines show little change above this 
column density, but the optical Fe 11 increases a factor 
^ 0.5 dex from ~ 10^^ cm~^ to ^ 10^^ cm~^. Using 
arguments based on viria l dete rminations of the black 
hole mass in AGN. lNetzeii ()2008D also concludes that the 
column densities must substantially exceed ~ 10^^ cm~^ 
to avoid excessive effects of radiation pressure on the 
orbital velocities of the BLR clouds. Thus, there may be 
limited freedom to vary the column density in order to 
produce the wide range of optical Fe 11 strength observed. 

The relative behavior of the optical and UV Fe 11 bands 
may provide clues to the predominant cause of the range 
of Fe II strength. In our models, the UV and optical Fe 11 
both increase with Fe abundance, albeit more weakly 
for the UV blends. However, increasing the microtur- 
bulence increased the UV Fe 11 by a greater factor that 
the optical Fe ii; and increasing the column density be- 
yond ~ 10^'^ cm~^ mainly i ncrea ses the optical Fe 11. as 
noted above. iShang et al.l (|2007[ ) give optical and Fe 11 
strengths for a sample of AGN. Their results show a much 
greater range in the EW of the optical Fc 11 bands than in 
the UV bands, and a wea k anti-correlation between th e 
optical and UV (see also iWills. Netzer. fc WiH^ [19851 ) . 
The fact that the optical Fe 11 shows a greater range 
of intensity that the UV Fe 11 may favor an explana- 
tion other than microturbulence. One complication is 
whether the optical and UV lines originate at substan- 
tially different radii in the BLR, as sugges ted by some 
reverberation a n d line- width s t udies (e.g., iMaod 119931 : 
IHu et al.|[2008l ). IFerland et all (|2009f) suggest that the 
observed optical Fe 11 may be strongly affected by radia- 
tion escaping from the shielded face of the photoionizcd 
clouds. The observational and theoretical situation for 
Fe II is complex, and further work will be needed to de- 
vise definitive tests of the role of chemical abundances 
and physical conditions. 

6. CONCLUSION 
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We have used composite SDSS spectra of AGN to ex- 
amine the behavior of the narrow emission hnes as a func- 
tion of Fe II strength. The [Fe Vii] Une shows only a weak 
increase with increasing Fc ii strength, indicating that 
the iron abundance contributes httle to the wide range 
of Fe II strength in AGN. There is, however, a significant 
increase in the N/0 abundance ratio with Fe ii strength, 
which suggests an increase in overall metallicity. There 
is little change in the level of ionization in the NLR as 
a function of Fc ii strength. This, together with results 
of photoionization models, suggests that differences in 
the shape of the ionizing continuum, specifically the soft 
X-ray luminosity, are not the main drivers of the Fe ii 
strength. We propose that differences in the degree of 
depletion of Fe into grains in the low ionization portion 
of the BLR are largely responsible for the weakness of 
Fe II in some AGN, while it is strong in others. This 
picture is consistent with the approximate coincidence of 
the BLR radius and the dust sublimation radius, with 
indications that the BLR consists of material flowing in- 
ward from the dusty torus toward the central black hole, 
and with the variation of grain temperature with depth 
in the ionization structure of the BLR gas. 

The strength of Fc ii emission is a major component of 
the set of correlations known as "Eigenvector 1" (EVl) 
discussed by BG92. These include weak [O iii] associ- 
ated with strong Fe ii and narrower widths of the broad 
H/3 line. Radio loud AGN tend to have strong [O ill] and 
weak Fc II. These trends have been the subject of many 
studies, but a good physical understanding of their origin 
remains lacking. BG92 suggested that high column den- 
sities in the BLR enhance Fe ii whi le diminishing t he ion - 
izing radiation reaching the NLR. iLudwig et al.l ()2009[ ). 
in a spectral principal components analysis of AGN in 
SDSS, found the interpretation of the eigenvectors to be 
complicated. They argued that covering factor of the 
NLR was the l ikely c ause of the range in [O ill] strength. 
iFerland et al.l ()2009D suggest that the higher column den- 
sities required for infall in more luminous AGN can help 
to explain the correlation of Fe ii strength with L/ L^^. 

The interpretation of Fe ii strength in terms of dust 



depletion opens many questions for investigation. Can 
support for this picture be found in the line intensities 
of other refractory elements? Does the Fe ii emitting 
radius bear a different relationship to the dust sublima- 
tion radius for strong and weak Fe ii emitters? What are 
the implications for the infrared emission of AGN? What 
underlying causes lead to the correlations between Fe ii 
strength and other properties such as [O ill] strength, ra- 
dio emission, and Eddington ratio? This paper does not 
offer answers to these larger questions, but the explana- 
tion of Fe II strength in terms of gas-phase depletions 
gives a new context in which to address them. 
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Fig. 1. — Composite SDSS spectra for the five bins in optical Fe II strength. Vertical axis gives specific flux Fx ■ Note the broad Fe II 
blends at A4570 and A5250. See text for discussion. 




vil] and [Nc v] narrow emission lines in the composite spectra. See text for discussion. 
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Fig. 3.— 



Narrow emission-line intensity ratios ratios for the five composite spectra. See text for discussion including errors. 



